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Description 

[TRANSIENT SIMULATION USING 
ADAPTIVE PIECEWISE CONSTANT 

MODEL] 

Background of Invention 
[0001] Technical Field 

[0002] The present invention relates generally to transient simu- 
lators, and more particularly, to a transient simulation 
system, methods and program product that implement an 
adaptive piecewise constant model. 

[0003] Related Art 

[0004] with the continuing miniaturization of very large scale in- 
tegrated (VLSI) circuitry, interest in conducting formal cir- 
cuit tuning using circuit transient simulation for gradient 
calculation is increasing. "Transient simulation" is the 
analysis of a circuit during time-varying conditions. "Tun- 
ing" includes optimization of a number of circuit parame- 
ters, e.g., device widths, to maximize circuit performance 



relative to circuit parameters. A challenge for circuit tran- 
sient simulation is improving simulator performance, i.e., 
run time. Tuning also requires conducting another simu- 
lation to determine "gradients" for the circuit, which are 
circuit measurement sensitivities to other circuit parame- 
ters, e.g., a transistor's width effect on circuit delay, resis- 
tor size effect on circuit speed, etc. This technique further 
increases the importance of performance. Performance of 
a simulator is based on how many operations are required 
to determine circuit performance and gradients. Unfortu- 
nately, it is common for a particular simulation operation 
to be called millions of times during a tuning run. 
[0005] one common simulation approach tests a circuit design 
by conducting simulations based on directly solving for 
voltages and currents via a simultaneous solve technique. 
Many of these type transient simulators provide for opti- 
mization of simulation by allowing an adaptive timing 
step, with the allowable timing step being chosen to en- 
sure that some user-defined error criteria is satisfied. For 
example, a local truncation error (LTE) may be imple- 
mented that describes a single-step voltage error as a 
function of the time step. This adaptive timing step is ad- 
vantageous for two reasons. First, it allows the simulator 



to adapt its rate of progress as needed to maintain accu- 
racy, while maximizing performance. Second, it allows the 
user to have some level of control over a performance/accu- 
racy tradeoff. There are some types of event-based fast 
simulators, however, for which this adaptive timing step 
capability is not possible. 
[0006] one type of simulator in which an adaptive timing step is 
not possible are those that rely on piecewise constant 
(PWC) table models to represent the simulation parame- 
ters. FIG. 1 illustrates a graphical representation of a PWC 
table model in which a change in current (AI r ) is constant 
within the model and a change in voltage (AV) is constant 
for each step within the PWC model. In some cases, a 
change in current (Al ) may vary. PWC table models know 

R 

current and voltage changes occurring at a particular cir- 
cuit element at issue, but do not know the timing of the 
changes. PWC simulators solve for the time to the next 
step (h), or "segment boundary," in the PWC model. Once 
the time solving step is complete, the relevant values are 
then obtained directly from a PWC table model, i.e., these 
parameters remain constant between steps by definition. 
This model is in contrast to the above-described simula- 
tors in which a timing change is known and the current 



and voltage parameter values are simultaneously deter- 
mined. 

[0007] Although PWC approximation aids in rapid circuit simula- 
tion by limiting all circuit changes to predefined values at 
discrete points in time, this approach also makes it diffi- 
cult to balance performance and accuracy available in 
other simulators because the integration scheme and sim- 
ulation accuracy are essentially predefined in the PWC ta- 
ble models. That is, the simulation steps are determined 
solely by the table step size. To ensure accurate simula- 
tion with these simulators, the step sizes of the PWC 
model must be made fine enough to maintain accuracy 
over the worst case predicted simulation conditions. Un- 
fortunately, this situation typically results in a very con- 
servative level of accuracy for at least a portion of the 
simulation so as not to lose performance. 

[0008] while it is possible to affect performance versus accuracy 
at a crude level in a PWC simulator, for example, through 
definition of a single optimized table step size during PWC 
model generation, it is not possible to implement, for ex- 
ample, an LTE relation to provide a controlled perfor- 
mance/accuracy balance. This inefficiency is particularly 
troublesome given the typical application for a PWC simu- 



lation approach - circuit tuning. Despite PWC simulators 
not providing the most efficient approach for transient 
simulation, their use remains advantageous because they 
provide the most efficient means of evaluating sensitivi- 
ties for non-stiff systems. Unfortunately, as gradients are 
typically calculated many times inside a particular tuning 
loop, the effect on performance of a PWC simulation is 
significant. 

[0009] | n V j ew 0 f the foregoing, there is a need in the art for way 

of improving simulation performance for a PWC-based 

transient simulator. 
Summary of Invention 

[0010] jhe invention includes a transient simulation system, 

methods and program product that implement an adap- 
tive piecewise constant (PWC) model. The invention evalu- 
ates an error criteria to determine a maximum allowable 
change in one of a current and a voltage; and simulates 
the transient conditions by implementing an adaptive step 
in the PWC model according to the maximum allowable 
change. The invention allows dynamic or static adaptation 
of a PWC model according to an error criteria. 

[001 1] a first aspect of the invention is directed to a method for 
simulating transient conditions in a circuit using a piece- 



wise constant model, the method comprising the steps of: 
evaluating an error criteria to determine a maximum al- 
lowable change in a current-voltage pair; and simulating 
the transient conditions by implementing an adaptive step 
in the piecewise constant model according to the maxi- 
mum allowable change. 
[0012] a second aspect of the invention is directed to a method 
for simulating transient conditions in a circuit using a 
piecewise constant model including a plurality of steps, 
the method comprising the steps of: replacing a plurality 
of first steps in the piecewise constant model with a lesser 
number of second steps to address an error criteria; and 
simulating the transient conditions using the piecewise 
constant model including the lesser number of second 
steps. 

[0013] a third aspect of the invention is directed to a computer 
program product comprising a computer useable medium 
having computer readable program code embodied 
therein for simulating transients conditions in a circuit us- 
ing a piecewise constant model, the program product 
comprising: program code configured to evaluate an error 
criteria to determine a maximum allowable change in a 
current-voltage pair; and program code configured to 



simulate the transient conditions by implementing an 
adaptive step in the piecewise constant model according 
to the maximum allowable change. 

[0014] a fourth aspect of the invention is directed to a system for 
simulating transient conditions in a circuit using a piece- 
wise constant model, the system comprising: means for 
evaluating an error criteria to determine a maximum al- 
lowable change in a current-voltage pair; and means for 
simulating the transient conditions by implementing an 
adaptive step in the piecewise constant model according 
to the maximum allowable change. 

[0015] The foregoing and other features of the invention will be 

apparent from the following more particular description of 

embodiments of the invention. 
Brief Description of Drawings 

[0016] The embodiments of this invention will be described in 

detail, with reference to the following figures, wherein like 
designations denote like elements, and wherein: 

[0017] FIG. 1 shows a graphical representation of a conventional 
piecewise constant (PWC) model for use by a simulation 
system. 

[0018] FIG. 2 shows a block diagram of a simulation system ac- 
cording to the invention. 



[0019] FIG. 3 shows a flow diagram of operation methodology 
according to a first embodiment of the invention. 

[0020] FIG. 4 shows a graphical representation of an adaptive 
PWC model. 

[0021] FIG. 5 shows a graphical representation of rendering of an 
adaptive step of an adaptive PWC model according to an 
optional embodiment of the invention. 

[0022] FIG. 6 shows a graphical representation comparing a con- 
ventional PWC model to an adaptive PWC model. 
Detailed Description 

[0023] For purposes of clarity only, the detailed description in- 
cludes the following outline: I. PWC-based Simulation Sys- 
tem Overview; II. Operational Methodology; and III. Con- 
clusion. 

[0024] The present invention includes a system, methods and 
program product to achieve a performance/accuracy 
tradeoff while using a piecewise constant (PWC) simulation 
approach. As the integration scheme of PWC-based simu- 
lators is essentially locked into the corresponding PWC 
models (as described above), the invention provides away 
to define adaptive PWC models, with the step size adapt- 
ing as needed to satisfy error criteria. The invention al- 
lows the simulation software to make small steps in criti- 



cal regions to maintain accuracy, while allowing larger 
steps where possible for performance. The invention also 
allows a single model to be used to represent all linear el- 
ements (such as resistors in the case of electrical simula- 
tion), as the optimal step size automatically adapts with 
element value. These benefits both improve the computa- 
tional efficiency of the code and simplify the definition of 
element models. 
[0025] |, pwC-based Simulation System Overview 

[0026] with reference to the accompanying drawings, FIG. 2 is a 
block diagram of a piecewise constant (PWC) based simu- 
lation system 10 in accordance with the invention. Simu- 
lation system 10 includes a memory 12, a processing unit 
(PU) 14, input/output devices (I/O) 16 and a bus 18. A 
database 20 may also be provided for storage of data rel- 
ative to processing tasks. Memory 12 includes a program 
product 22 that, when executed by PU 14, comprises vari- 
ous functional capabilities described in further detail be- 
low. Memory 12 (and database 20) may comprise any 
known type of data storage system and/or transmission 
media, including magnetic media, optical media, random 
access memory (RAM), read only memory (ROM), a data 
object, etc. Moreover, memory 12 (and database 20) may 



reside at a single physical location comprising one or 
more types of data storage, or be distributed across a 
plurality of physical systems. PU 14 may likewise comprise 
a single processing unit, or a plurality of processing units 
distributed across one or more locations. I/O 16 may 
comprise any known type of input/output device including 
a network system, modem, keyboard, mouse, scanner, 
voice recognition system, CRT, printer, disc drives, etc. 
Additional components, such as cache memory, commu- 
nication systems, system software, etc., may also be in- 
corporated into simulation system 10. 

[0027] As shown in FIG. 2, program product 22 may include an 
evaluator 24 including a change calculator 26, a replacer 
28, an adaptive step limiter 30, and an adaptive step Ten- 
derer 32; a simulator 40 and other system components 
50. Simulator 40 includes any now known or later devel- 
oped PWC-based simulator. Other system components 50 
include any other function necessary or desired for simu- 
lation system 10 not otherwise explicitly described herein. 

[0028] As also shown in FIG. 2, an error criteria 62 may be input 
to simulation system 10, as will be described in more de- 
tail below. Database 20 includes a number of PWC models 
60 that may be recalled and then modified as described 



herein below. 
[0029] ||. Operational Methodology 

[0030] Turning to FIG. 3, a method for simulating transient con- 
ditions in a circuit using an adaptive PWC model is illus- 
trated according to a first embodiment of the invention. 

[0031] | n a fj rs t step SI, an error criteria 62 (FIG. 2) is evaluated 
to determine a maximum allowable change in a current- 
voltage pair in a PWC model 60 by evaluator 24 (FIG. 2). 
As used herein, "current-voltage pair" signifies a current- 
voltage relationship. A "change in a current-voltage pair" 
may be initiated in terms of a change in current or a 
change in voltage, but inherently results in a change in 
both. FIG. 4 illustrates a graphical representation of a PWC 
model 60 (FIG. 2) recalled from database 20. As shown, 
PWC model 60 includes equally proportioned steps 74 in 
which a current change (Al ) is constant and a voltage 

R 

change (AV) is constant for each step. PWC model 60 im- 
plements a plurality of steps representing integration step 
sizes. For a given voltage range, a current is maintained 
constant such that a graph of PWC model 60 appears in a 
stair form for models with liner current-voltage relations. 
For non-linear models, a voltage change remains constant 
while a current change adjusts as needed to track the 



function. Performance of PWC model 60, i.e., run time, is 
based on the number of steps that must be performed. 
Accordingly, fewer steps results in greater performance, 
but less accuracy. 
[0032] Overlayed over PWC model 60 is an adaptive PWC model 

70, according to the invention, in which an adaptation in a 
current-voltage pair is made (i.e, an adaptation of a cur- 
rent change (Al ) or a voltage change (AV)) within the 

R 

model to address error criteria 62. As illustrated, a plural- 
ity of predefined (first) steps 74 of PWC model 60 have 
been replaced by replacer 28 (FIG. 2) by a lesser number 
of (second) steps 76 of adaptive PWC model 70. Accord- 
ingly, adaptive PWC model 70 increases the size of steps, 
which reduces event count and reduces simulation time. 
An adaptive PWC model 70 may be implemented for each 
element of a circuit to be simulated. Adaptive PWC model 
70 can be implemented during runtime as needed in a dy- 
namic approach, or as a static approach, as will be de- 
scribed further below. 
[0033] Returning to FIG. 3 in conjunction with FIG. 2, error crite- 
ria 62 may be set to address a number of error-based sit- 
uations. In one preferred embodiment, error criteria 62 is 
based on a variety of an approximate relative timing error 



(At/t). An error criteria 62 can be input in a variety of 
forms. In one embodiment, an error criteria 62 may be in- 
put as a percentage, e.g., 5%. In this case, if an input to a 
node is expected to take 100 picoseconds (ps), and takes 
94 ps, then an error criteria of 5% would be violated, 
which would allow implementation of an adaptive step 76 
(FIG. 4), e.g., lengthening of a step by changing AV. An 
approximate relative timing error (At/t) can be defined as: 
[0034] At/t = (A |/v' c ) 2 

eq 

[0035] where At is time change between an initiating event (e.g., 
v(t) in FIG. 4 and a subsequent event v(t+h), i.e., At = 
(t+h) - t), t is time of the initiating event v(t), V is the 
derivative voltage for the present element, e.g., how fast 
voltage rising/falling across resistor (v/s)), Al is the cur- 
rent change, and C is an equivalent series capacitance 

eq 

across an element, which is most likely caused by device 
parasitics. 

[0036] a maximum allowable change in terms of current Al is 

max 

determined by change calculator 26 (FIG. 2) according to 

the rearrangement of the above equation to: 

[0037] Al = (At/t) 1/2 V C 

max eq 

[0038] while, in this case, the maximum allowable change is de- 



termined in terms of a current, it should be recognized 
that a maximum allowable change in terms of voltage (AV) 
can be derived from the current change (Al) using the 
known current-voltage pair relationship in a trapezoidal 
integration scheme (FIG. 5). This is represented as a "fine- 
grained" PWC nominal step-size model at the minimum 
allowable AV step, multiples of which may actually be 
taken based on Al (multiple = N , as described be- 

max step 

low). 

[0039] As shown in FIG. 3, a second step S2 includes simulating 
transient conditions by simulator 40 (FIG. 2) by imple- 
menting an adaptive step in the piecewise constant model 
according to the maximum allowable change. As noted 
above, simulator 40 (FIG. 2) may include any now known 
or later developed PWC-based simulator. 

[0040] with further regard to evaluation step SI, one challenge 
for the above evaluation steps is that simulation system 
10 evaluates each element separately. In this case, it is 
possible for an element to cause a large enough time 
change At in an adaptive step such that evaluation of 
other effects from neighboring elements would be over- 
looked. That is, a time of evaluation would skip necessary 
evaluations. In order to address this situation, as an op- 



tional sub-step S1A of evaluation step SI, the number of 
adaptive steps that adaptive PWC model 70 (FIG. 3) can 
include may be limited by adaptive step limiter 30 (FIG. 2). 
Sub-step S1A may be initiated by determining a number 
of steps N by dividing the maximum allowable change 
by a corresponding predefined minimum allowable 
change. That is, according to: 

[0041] N =A | /Al 

step max nominal 

[0042] where Al is the maximum allowable change in current 

max 

and Al is a nominal resolution, i.e, minimum allow- 

nominal 

able change. The number of steps N may also be cal- 
culated in terms of voltage. Calculation of N could re- 
step 

12 

suit in a number ranging from 0 to le . Where N is a 

step 

very large number (e.g., > 3000), implementation may not 
be advisable. To prevent overburdening simulation system 
10, the number of steps N may be limited to some "rea- 

step 

sonable" value by adaptive step limiter 30 (FIG. 2). For ex- 
ample, the following logic may be implemented to limit 
the number of steps: if N < 1, then N =1; else, if N 

step step step 

< 100, then N = 100; else, if N < 3000, then N 

step step step 

= 500; else, N = 750. It should be recognized that dif— 

step 

ferent situations may require different ranges of the num- 
ber of steps, and that the above description is only i 1 1 us— 



trative. 

[0043] | n another optional sub-step SIB, the adaptive step may 
be rejected by replacer 28 (FIG. 2) when the sign of the 
derivative voltage across a circuit element of interest re- 
verses. Conventional PWC simulators never allow an in- 
stantaneous change in sign on a derivative voltage as it 
can result in numerical instability. Simulation system 10 
provides the same safety feature by rejecting the value of 
N calculated by the adaptive feature and forcing N 

step step 

to equal 1 if the derivative voltage reverses sign, e.g., 
positive to negative. 
[0044] | n another optional sub-step SIC, an adaptive step is ren- 
dered, i.e., inserted into PWC model 70, by adaptive step 
renderer 32 (FIG. 2) at an average value of the maximum 
allowable change. In particular, as shown in FIG. 5, if a 
maximum allowable change Al provides any value 

max 

other than a nominal Al (minimum), the adaptive step 174 
must not be exhausted on the current step. In particular, 
it is most advantageous if the adaptive step crosses a 
middle of an original step to maintain a trapezoidal inte- 
gration scheme. As shown in FIG. 5, adaptive step 174 is 
implemented over two subsequent steps to maintain a 
trapezoidal scheme. 



[0045] Dynamic and Static Models 

[0046] | n one embodiment, the above-described evaluating step 
SI executes dynamically during the simulating step S2 
such that for each element at least one adaptive step can 
be implemented, i.e, adaptive PWC model 70 (FIG. 3) is 
continually changing. FIG. 6 illustrates a graphical repre- 
sentation comparing a conventional PWC model to an 
adaptive PWC model according to the dynamic approach. 
This dynamic approach allows for the most adaptability in 
the PWC model definition, as each step is defined to sat- 
isfy a current simulation condition. However, to leverage 
this adaptability, a very fine nominal step size is required. 
Also, as the adaptive algorithm must be executed at every 
simulation step, it also can carry a moderately high com- 
putational overhead. This is particularly true for nonlinear 
element models, where it is not possible to interpolate the 
element parameters over multiple steps. Therefore, this 
approach is generally better suited for linear elements. 
This linearity allows the element "model" to consist of a 
single segment definition, and as the model consists of 
single real-valued piece of data, there is no overhead to 
making this segment relatively small (as opposed to mod- 
els in which the memory requirements go up with the di- 



mensionality of the element parameters). The actual step 
size required during any simulation step is then calculated 
to satisfy the error requirement, and the independent 
variable is linearly interpolated from this step size. 

[0047] However, in an alternative embodiment, the evaluating 

step SI executes prior to the simulating step S2 such that 
adaptive PWC model 70 (FIG. 3) is static. In the static ap- 
proach, the adaptive PWC model is predefined, and re- 
mains constant during simulation. While this approach has 
the drawback of not being able to react to the changing 
conditions that exist during simulation, it also incurs no 
computational overhead during simulation. The static ap- 
proach is also less dependent on the nominal segment 
size used in the original PWC model. Therefore, the static 
approach can be slightly more conservative, but is gener- 
ally applicable to all circuit elements. 

[0048] III. Conclusion 

[0049] | n the previous discussion, it will be understood that the 
method steps discussed are performed by a processor, 
such as PU 14 of system 10, executing instructions of 
program product 22 stored in memory. It is understood 
that the various devices, modules, mechanisms and sys- 
tems described herein may be realized in hardware, soft- 



ware, or a combination of hardware and software, and 
may be compartmentalized other than as shown. They 
may be implemented by any type of computer system or 
other apparatus adapted for carrying out the methods de- 
scribed herein. A typical combination of hardware and 
software could be a general-purpose computer system 
with a computer program that, when loaded and exe- 
cuted, controls the computer system such that it carries 
out the methods described herein. Alternatively, a specific 
use computer, containing specialized hardware for carry- 
ing out one or more of the functional tasks of the inven- 
tion, could be utilized. The present invention can also be 
embedded in a computer program product, which com- 
prises all the features enabling the implementation of the 
methods and functions described herein, and which - 
when loaded in a computer system - is able to carry out 
these methods and functions. Computer program, soft- 
ware program, program, program product, or software, in 
the present context mean any expression, in any lan- 
guage, code or notation, of a set of instructions intended 
to cause a system having an information processing capa- 
bility to perform a particular function either directly or af- 
ter the following: (a) conversion to another language, code 



or notation; and/or (b) reproduction in a different material 
form. 

[0050] while this invention has been described in conjunction 
with the specific embodiments outlined above, it is evi- 
dent that many alternatives, modifications and variations 
will be apparent to those skilled in the art. Accordingly, 
the embodiments of the invention as set forth above are 
intended to be illustrative, not limiting. Various changes 
may be made without departing from the spirit and scope 
of the invention as defined in the following claims. 



